AmeriFlux BASE data pipeline to support network growth and data sharing.

AmeriFlux is a network of research sites that measure carbon, water, and energy fluxes between ecosystems and the atmosphere using the eddy covariance technique to study a variety of Earth science questions. AmeriFlux's diversity of ecosystems, instruments, and data-processing routines create challenges for data standardization, quality assurance, and sharing across the network. To address these challenges, the AmeriFlux Management Project (AMP) designed and implemented the BASE data-processing pipeline. The pipeline begins with data uploaded by the site teams, followed by the AMP team's quality assurance and quality control (QA/QC), ingestion of site metadata, and publication of the BASE data product. The semi-automated pipeline enables us to keep pace with the rapid growth of the network. As of 2022, the AmeriFlux BASE data product contains 3,130 site years of data from 444 sites, with standardized units and variable names of more than 60 common variables, representing the largest long-term data repository for flux-met data in the world. The standardized, quality-ensured data product facilitates multisite comparisons, model evaluations, and data syntheses.

Detection and attribution of an anomaly in terrestrial photosynthesis in Europe during the COVID-19 lockdown.

Carbon dioxide (CO2) uptake by plant photosynthesis, referred to as gross primary production (GPP) at the ecosystem level, is sensitive to environmental factors, including pollutant exposure, pollutant uptake, and changes in the scattering of solar shortwave irradiance (SWin) - the energy source for photosynthesis. The 2020 spring lockdown due to COVID-19 resulted in improved air quality and atmospheric transparency, providing a unique opportunity to assess the impact of air pollutants on terrestrial ecosystem functioning. However, detecting these effects can be challenging as GPP is influenced by other meteorological drivers and management practices. Based on data collected from 44 European ecosystem-scale CO2 flux monitoring stations, we observed significant changes in spring GPP at 34 sites during 2020 compared to 2015-2019. Among these, 14 sites showed an increase in GPP associated with higher SWin, 10 sites had lower GPP linked to atmospheric and soil dryness, and seven sites were subjected to management practices. The remaining three sites exhibited varying dynamics, with one experiencing colder and rainier weather resulting in lower GPP, and two showing higher GPP associated with earlier spring melts. Analysis using the regional atmospheric chemical transport model (LOTOS-EUROS) indicated that the ozone (O3) concentration remained relatively unchanged at the research sites, making it unlikely that O3 exposure was the dominant factor driving the primary production anomaly. In contrast, SWin increased by 9.4 % at 36 sites, suggesting enhanced GPP possibly due to reduced aerosol optical depth and cloudiness. Our findings indicate that air pollution and cloudiness may weaken the terrestrial carbon sink by up to 16 %. Accurate and continuous ground-based observations are crucial for detecting and attributing subtle changes in terrestrial ecosystem functioning in response to environmental and anthropogenic drivers.

A widely-used eddy covariance gap-filling method creates systematic bias in carbon balance estimates.

Climate change mitigation requires, besides reductions in greenhouse gas emissions, actions to increase carbon sinks in terrestrial ecosystems. A key measurement method for quantifying such sinks and calibrating models is the eddy covariance technique, but it requires imputation, or gap-filling, of missing data for determination of annual carbon balances of ecosystems. Previous comparisons of gap-filling methods have concluded that commonly used methods, such as marginal distribution sampling (MDS), do not have a significant impact on the carbon balance estimate. By analyzing an extensive, global data set, we show that MDS causes significant carbon balance errors for northern (latitude [Formula: see text]) sites. MDS systematically overestimates the carbon dioxide (CO[Formula: see text]) emissions of carbon sources and underestimates the CO[Formula: see text] sequestration of carbon sinks. We also reveal reasons for these biases and show how a machine learning method called extreme gradient boosting or a modified implementation of MDS can be used to substantially reduce the northern site bias.

Intercomparison of methods to estimate gross primary production based on CO2 and COS flux measurements.

Separating the components of ecosystem-scale carbon exchange is crucial in order to develop better models and future predictions of the terrestrial carbon cycle. However, there are several uncertainties and unknowns related to current photosynthesis estimates. In this study, we evaluate four different methods for estimating photosynthesis at a boreal forest at the ecosystem scale, of which two are based on carbon dioxide (CO2) flux measurements and two on carbonyl sulfide (COS) flux measurements. The CO2-based methods use traditional flux partitioning and artificial neural networks to separate the net CO2 flux into respiration and photosynthesis. The COS-based methods make use of a unique 5-year COS flux data set and involve two different approaches to determine the leaf-scale relative uptake ratio of COS and CO2 (LRU), of which one (LRUCAP) was developed in this study. LRUCAP was based on a previously tested stomatal optimization theory (CAP), while LRUPAR was based on an empirical relation to measured radiation. For the measurement period 2013-2017, the artificial neural network method gave a GPP estimate very close to that of traditional flux partitioning at all timescales. On average, the COS-based methods gave higher GPP estimates than the CO2-based estimates on daily (23% and 7% higher, using LRUPAR and LRUCAP, respectively) and monthly scales (20% and 3% higher), as well as a higher cumulative sum over 3 months in all years (on average 25% and 3% higher). LRUCAP was higher than LRU estimated from chamber measurements at high radiation, leading to underestimation of midday GPP relative to other GPP methods. In general, however, use of LRUCAP gave closer agreement with CO2-based estimates of GPP than use of LRUPAR. When extended to other sites, LRUCAP may be more robust than LRUPAR because it is based on a physiological model whose parameters can be estimated from simple measurements or obtained from the literature. In contrast, the empirical radiation relation in LRUPAR may be more site-specific. However, this requires further testing at other measurement sites.

Global nature run data with realistic high-resolution carbon weather for the year of the Paris Agreement.

The CO2 Human Emissions project has generated realistic high-resolution 9 km global simulations for atmospheric carbon tracers referred to as nature runs to foster carbon-cycle research applications with current and planned satellite missions, as well as the surge of in situ observations. Realistic atmospheric CO2, CH4 and CO fields can provide a reference for assessing the impact of proposed designs of new satellites and in situ networks and to study atmospheric variability of the tracers modulated by the weather. The simulations spanning 2015 are based on the Copernicus Atmosphere Monitoring Service forecasts at the European Centre for Medium Range Weather Forecasts, with improvements in various model components and input data such as anthropogenic emissions, in preparation of a CO2 Monitoring and Verification Support system. The relative contribution of different emissions and natural fluxes towards observed atmospheric variability is diagnosed by additional tagged tracers in the simulations. The evaluation of such high-resolution model simulations can be used to identify model deficiencies and guide further model improvements.

Direct observations of CO2 emission reductions due to COVID-19 lockdown across European urban districts.

The measures taken to contain the spread of COVID-19 in 2020 included restrictions of people's mobility and reductions in economic activities. These drastic changes in daily life, enforced through national lockdowns, led to abrupt reductions of anthropogenic CO2 emissions in urbanized areas all over the world. To examine the effect of social restrictions on local emissions of CO2, we analysed district level CO2 fluxes measured by the eddy-covariance technique from 13 stations in 11 European cities. The data span several years before the pandemic until October 2020 (six months after the pandemic began in Europe). All sites showed a reduction in CO2 emissions during the national lockdowns. The magnitude of these reductions varies in time and space, from city to city as well as between different areas of the same city. We found that, during the first lockdowns, urban CO2 emissions were cut with respect to the same period in previous years by 5% to 87% across the analysed districts, mainly as a result of limitations on mobility. However, as the restrictions were lifted in the following months, emissions quickly rebounded to their pre-COVID levels in the majority of sites.

Substantial hysteresis in emergent temperature sensitivity of global wetland CH4 emissions.

Wetland methane (CH4) emissions ([Formula: see text]) are important in global carbon budgets and climate change assessments. Currently, [Formula: see text] projections rely on prescribed static temperature sensitivity that varies among biogeochemical models. Meta-analyses have proposed a consistent [Formula: see text] temperature dependence across spatial scales for use in models; however, site-level studies demonstrate that [Formula: see text] are often controlled by factors beyond temperature. Here, we evaluate the relationship between [Formula: see text] and temperature using observations from the FLUXNET-CH4 database. Measurements collected across the globe show substantial seasonal hysteresis between [Formula: see text] and temperature, suggesting larger [Formula: see text] sensitivity to temperature later in the frost-free season (about 77% of site-years). Results derived from a machine-learning model and several regression models highlight the importance of representing the large spatial and temporal variability within site-years and ecosystem types. Mechanistic advancements in biogeochemical model parameterization and detailed measurements in factors modulating CH4 production are thus needed to improve global CH4 budget assessments.

Modeling Ambitions Outpace Observations of Forest Carbon Allocation.

There have been vociferous calls for 'tree-centered' vegetation models to refine predictions of forest carbon (C) cycling. Unfortunately, our global survey at flux-tower sites indicates insufficient empirical data support for this much-needed model development. We urge for a new generation of studies across large environmental gradients that strategically pair long-term ecosystem monitoring with manipulative experiments on mature trees. For this, we outline a versatile experimental framework to build cross-scale data archives of C uptake and allocation to structural, non-structural, and respiratory sinks. Community-wide efforts and discussions are needed to implement this framework, especially in hitherto underrepresented tropical forests. Global coordination and realistic priorities for data collection will thereby be key to achieve and maintain adequate empirical support for tree-centered vegetation modeling.

Partitioning net carbon dioxide fluxes into photosynthesis and respiration using neural networks.

The eddy covariance (EC) technique is used to measure the net ecosystem exchange (NEE) of CO2 between ecosystems and the atmosphere, offering a unique opportunity to study ecosystem responses to climate change. NEE is the difference between the total CO2 release due to all respiration processes (RECO), and the gross carbon uptake by photosynthesis (GPP). These two gross CO2 fluxes are derived from EC measurements by applying partitioning methods that rely on physiologically based functional relationships with a limited number of environmental drivers. However, the partitioning methods applied in the global FLUXNET network of EC observations do not account for the multiple co-acting factors that modulate GPP and RECO flux dynamics. To overcome this limitation, we developed a hybrid data-driven approach based on combined neural networks (NNC-part ). NNC-part incorporates process knowledge by introducing a photosynthetic response based on the light-use efficiency (LUE) concept, and uses a comprehensive dataset of soil and micrometeorological variables as fluxes drivers. We applied the method to 36 sites from the FLUXNET2015 dataset and found a high consistency in the results with those derived from other standard partitioning methods for both GPP (R2  > .94) and RECO (R2  > .8). High consistency was also found for (a) the diurnal and seasonal patterns of fluxes and (b) the ecosystem functional responses. NNC-part performed more realistic than the traditional methods for predicting additional patterns of gross CO2 fluxes, such as: (a) the GPP response to VPD, (b) direct effects of air temperature on GPP dynamics, (c) hysteresis in the diel cycle of gross CO2 fluxes, (d) the sensitivity of LUE to the diffuse to direct radiation ratio, and (e) the post rain respiration pulse after a long dry period. In conclusion, NNC-part is a valid data-driven approach to provide GPP and RECO estimates and complementary to the existing partitioning methods.

Monitoring tropical forests under a functional perspective with satellite-based vegetation optical depth.

Monitoring ecosystem functions in forests is a priority in a climate change scenario, as climate-induced events may initially alter the functions more than slow-changing attributes, such as biomass. The ecosystem functional properties (EFPs) are quantities that characterize key ecosystem processes. They can be derived by point observations of gas and energy exchanges between the ecosystems and the atmosphere that are collected globally at FLUXNET flux tower sites and upscaled at ecosystem level. The properties here considered describe the ability of ecosystems to optimize the use of resources for carbon uptake. They represent functional forest information, are dependent on environmental drivers, linked to leaf traits and forest structure, and influenced by climate change effects. The ability of vegetation optical depth (VOD) to provide forest functional information is investigated using 2011-2014 satellite data collected by the Soil Moisture and Ocean Salinity mission and using the EFPs as reference dataset. Tropical forests in Africa and South America were analyzed, also according to ecological homogeneous units. VOD jointly with water deficit information explained 93% and 87% of the yearly variability in both flux upscaled maximum gross primary productivity and light use efficiency functional properties, in Africa and South America forests respectively. Maps of the retrieved properties evidenced changes in forest functional responses linked to anomalous climate-induced events during the study period. The findings indicate that VOD can support the flux upscaling process in the tropical range, affected by high uncertainty, and the detection of forest anomalous functional responses. Preliminary temporal analysis of VOD and EFP signals showed fine-grained variability in periodicity, in signal dephasing, and in the strength of the relationships. In selected drier forest types, these satellite data could also support the monitoring of functional dynamics.

The FLUXCOM ensemble of global land-atmosphere energy fluxes.

Although a key driver of Earth's climate system, global land-atmosphere energy fluxes are poorly constrained. Here we use machine learning to merge energy flux measurements from FLUXNET eddy covariance towers with remote sensing and meteorological data to estimate global gridded net radiation, latent and sensible heat and their uncertainties. The resulting FLUXCOM database comprises 147 products in two setups: (1) 0.0833° resolution using MODIS remote sensing data (RS) and (2) 0.5° resolution using remote sensing and meteorological data (RS + METEO). Within each setup we use a full factorial design across machine learning methods, forcing datasets and energy balance closure corrections. For RS and RS + METEO setups respectively, we estimate 2001-2013 global (±1 s.d.) net radiation as 75.49 ± 1.39 W m-2 and 77.52 ± 2.43 W m-2, sensible heat as 32.39 ± 4.17 W m-2 and 35.58 ± 4.75 W m-2, and latent heat flux as 39.14 ± 6.60 W m-2 and 39.49 ± 4.51 W m-2 (as evapotranspiration, 75.6 ± 9.8 × 103 km3 yr-1 and 76 ± 6.8 × 103 km3 yr-1). FLUXCOM products are suitable to quantify global land-atmosphere interactions and benchmark land surface model simulations.

The ecosystem carbon sink implications of mountain forest expansion into abandoned grazing land: The role of subsoil and climatic factors.

Woody encroachment is a widespread phenomenon resulting from the abandonment of mountain agricultural and pastoral practices during the last century. As a result, forests have expanded, increasing biomass and necromass carbon (C) pools. However, the impact on soil organic carbon (SOC) is less clear. The main aim of this study was to investigate the effect of woody encroachment on SOC stocks and ecosystem C pools in six chronosequences located along the Italian peninsula, three in the Alps and three in the Apennines. Five stages along the chronosequences were identified in each site. Considering the topsoil (0-30 cm), subsoil (30 cm-bedrock) and whole soil profile, the temporal trend in SOC stocks was similar in all sites, with an initial increment and subsequent decrement in the intermediate phase. However, the final phase of the woody encroachment differed significantly between the Alps (mainly conifers) and the Apennines (broadleaf forests) sites, with a much more pronounced increment in the latter case. Compared to the previous pastures, after mature forest (>62 years old) establishment, SOC stocks increased by: 2.1(mean) ±â€¯18.1(sd) and 50.1 ±â€¯25.2 Mg C·ha-1 in the topsoil, 7.3 ±â€¯17.4 and 93.2 ±â€¯29.7 Mg C·ha-1 in the subsoil, and 9.4 ±â€¯24.4 and 143.3 ±â€¯51.0 Mg C·ha-1 in the whole soil profile in Alps and Apennines, respectively. Changes in SOC stocks increased with mean annual air temperature and average minimum winter temperature, and were negatively correlated with the sum of summer precipitation. At the same time, all other C pools (biomass and necromass) increased by 179.1 ±â€¯51.3 and 304.2 ±â€¯67.6 Mg C·ha-1 in the Alps and the Apennines sites, respectively. This study highlights the importance of considering both the subsoil, since deep soil layers contributed 38% to the observed variations in carbon stocks after land use change, and the possible repercussions for the carbon balance of large areas where forests are expanding, especially under pressing global warming scenarios.

Widespread inhibition of daytime ecosystem respiration.

The global land surface absorbs about a third of anthropogenic emissions each year, due to the difference between two key processes: ecosystem photosynthesis and respiration. Despite the importance of these two processes, it is not possible to measure either at the ecosystem scale during the daytime. Eddy-covariance measurements are widely used as the closest 'quasi-direct' ecosystem-scale observation from which to estimate ecosystem photosynthesis and respiration. Recent research, however, suggests that current estimates may be biased by up to 25%, due to a previously unaccounted for process: the inhibition of leaf respiration in the light. Yet the extent of inhibition remains debated, and implications for estimates of ecosystem-scale respiration and photosynthesis remain unquantified. Here, we quantify an apparent inhibition of daytime ecosystem respiration across the global FLUXNET eddy-covariance network and identify a pervasive influence that varies by season and ecosystem type. We develop partitioning methods that can detect an apparent ecosystem-scale inhibition of daytime respiration and find that diurnal patterns of ecosystem respiration might be markedly different than previously thought. The results call for the re-evaluation of global terrestrial carbon cycle models and also suggest that current global estimates of photosynthesis and respiration may be biased, some on the order of magnitude of anthropogenic fossil fuel emissions.

Compensatory water effects link yearly global land CO2 sink changes to temperature.

Large interannual variations in the measured growth rate of atmospheric carbon dioxide (CO2) originate primarily from fluctuations in carbon uptake by land ecosystems. It remains uncertain, however, to what extent temperature and water availability control the carbon balance of land ecosystems across spatial and temporal scales. Here we use empirical models based on eddy covariance data and process-based models to investigate the effect of changes in temperature and water availability on gross primary productivity (GPP), terrestrial ecosystem respiration (TER) and net ecosystem exchange (NEE) at local and global scales. We find that water availability is the dominant driver of the local interannual variability in GPP and TER. To a lesser extent this is true also for NEE at the local scale, but when integrated globally, temporal NEE variability is mostly driven by temperature fluctuations. We suggest that this apparent paradox can be explained by two compensatory water effects. Temporal water-driven GPP and TER variations compensate locally, dampening water-driven NEE variability. Spatial water availability anomalies also compensate, leaving a dominant temperature signal in the year-to-year fluctuations of the land carbon sink. These findings help to reconcile seemingly contradictory reports regarding the importance of temperature and water in controlling the interannual variability of the terrestrial carbon balance. Our study indicates that spatial climate covariation drives the global carbon cycle response.

Above-ground woody carbon sequestration measured from tree rings is coherent with net ecosystem productivity at five eddy-covariance sites.

• Attempts to combine biometric and eddy-covariance (EC) quantifications of carbon allocation to different storage pools in forests have been inconsistent and variably successful in the past. • We assessed above-ground biomass changes at five long-term EC forest stations based on tree-ring width and wood density measurements, together with multiple allometric models. Measurements were validated with site-specific biomass estimates and compared with the sum of monthly CO2 fluxes between 1997 and 2009. • Biometric measurements and seasonal net ecosystem productivity (NEP) proved largely compatible and suggested that carbon sequestered between January and July is mainly used for volume increase, whereas that taken up between August and September supports a combination of cell wall thickening and storage. The inter-annual variability in above-ground woody carbon uptake was significantly linked with wood production at the sites, ranging between 110 and 370 g C m(-2) yr(-1) , thereby accounting for 10-25% of gross primary productivity (GPP), 15-32% of terrestrial ecosystem respiration (TER) and 25-80% of NEP. • The observed seasonal partitioning of carbon used to support different wood formation processes refines our knowledge on the dynamics and magnitude of carbon allocation in forests across the major European climatic zones. It may thus contribute, for example, to improved vegetation model parameterization and provides an enhanced framework to link tree-ring parameters with EC measurements.

Climate extremes and the carbon cycle.

The terrestrial biosphere is a key component of the global carbon cycle and its carbon balance is strongly influenced by climate. Continuing environmental changes are thought to increase global terrestrial carbon uptake. But evidence is mounting that climate extremes such as droughts or storms can lead to a decrease in regional ecosystem carbon stocks and therefore have the potential to negate an expected increase in terrestrial carbon uptake. Here we explore the mechanisms and impacts of climate extremes on the terrestrial carbon cycle, and propose a pathway to improve our understanding of present and future impacts of climate extremes on the terrestrial carbon budget.

Recent decline in the global land evapotranspiration trend due to limited moisture supply.

More than half of the solar energy absorbed by land surfaces is currently used to evaporate water. Climate change is expected to intensify the hydrological cycle and to alter evapotranspiration, with implications for ecosystem services and feedback to regional and global climate. Evapotranspiration changes may already be under way, but direct observational constraints are lacking at the global scale. Until such evidence is available, changes in the water cycle on land−a key diagnostic criterion of the effects of climate change and variability−remain uncertain. Here we provide a data-driven estimate of global land evapotranspiration from 1982 to 2008, compiled using a global monitoring network, meteorological and remote-sensing observations, and a machine-learning algorithm. In addition, we have assessed evapotranspiration variations over the same time period using an ensemble of process-based land-surface models. Our results suggest that global annual evapotranspiration increased on average by 7.1 ± 1.0 millimetres per year per decade from 1982 to 1997. After that, coincident with the last major El Niño event in 1998, the global evapotranspiration increase seems to have ceased until 2008. This change was driven primarily by moisture limitation in the Southern Hemisphere, particularly Africa and Australia. In these regions, microwave satellite observations indicate that soil moisture decreased from 1998 to 2008. Hence, increasing soil-moisture limitations on evapotranspiration largely explain the recent decline of the global land-evapotranspiration trend. Whether the changing behaviour of evapotranspiration is representative of natural climate variability or reflects a more permanent reorganization of the land water cycle is a key question for earth system science.

On the differential advantages of evergreenness and deciduousness in mediterranean oak woodlands: a flux perspective.

We assessed the differential advantages of deciduousness and evergreenness by examining 26 site-years of carbon dioxide, water vapor, and energy flux measurements from five comparable oak woodlands in France, Italy, Portugal, and California (USA). On average, the evergreen and deciduous oak woodlands assimilated and respired similar amounts of carbon while using similar amounts of water. These results suggest that evergreen and deciduous woodlands have specific, and similar, ecological costs in mediterranean climates, and that both leaf habits are able to meet these costs. What are the mechanisms behind these findings? Deciduous oaks compensated for having a shorter growing season by attaining a greater capacity to assimilate carbon for a given amount of intercepted solar radiation during the well-watered spring period; at saturating light levels, deciduous oaks gained carbon at six times the rate of evergreen oaks. Otherwise, the two leaf habits experienced similar efficiencies in carbon use (the change in carbon respired per change in carbon assimilated), water use (the change in carbon assimilation per change in water evaporated), and rainfall use (the change in evaporation per change in rainfall). Overall, leaf area index, rather than leaf habit, was the significant factor in determining the absolute magnitude of carbon gained and water lost by each evergreen and deciduous oak woodland over an annual interval; the closed canopies assimilated and respired more carbon and transpired more water than the open canopies. Both deciduous and evergreen mediterranean oaks survive in their seasonally hot/dry, wet/ cool native range by ensuring that actual evaporation is less than the supply of water. This feat is accomplished by adjusting the leaf area index to reduce total water loss at the landscape scale, by down-regulating photosynthesis, respiration, and stomatal conductance with progressive seasonal soil water deficits, and by extending their root systems to tap groundwater.

Terrestrial gross carbon dioxide uptake: global distribution and covariation with climate.

Terrestrial gross primary production (GPP) is the largest global CO(2) flux driving several ecosystem functions. We provide an observation-based estimate of this flux at 123 +/- 8 petagrams of carbon per year (Pg C year(-1)) using eddy covariance flux data and various diagnostic models. Tropical forests and savannahs account for 60%. GPP over 40% of the vegetated land is associated with precipitation. State-of-the-art process-oriented biosphere models used for climate predictions exhibit a large between-model variation of GPP's latitudinal patterns and show higher spatial correlations between GPP and precipitation, suggesting the existence of missing processes or feedback mechanisms which attenuate the vegetation response to climate. Our estimates of spatially distributed GPP and its covariation with climate can help improve coupled climate-carbon cycle process models.

ASPIS, A Flexible Multispectral System for Airborne Remote Sensing Environmental Applications.

Airborne multispectral and hyperspectral remote sensing is a powerful tool for environmental monitoring applications. In this paper we describe a new system (ASPIS) composed by a 4-CCD spectral sensor, a thermal IR camera and a laser altimeter that is mounted on a flexible Sky-Arrow airplane. A test application of the multispectral sensor to estimate durum wheat quality is also presented.